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A modified renormalization group equation for the inverse 
extrapolation length c is derived by considering the phase 
shifts of order parameter fluctuations. The resulting non- 
linear equation is also derived using standard methods and 
some additional assumptions. The associated renormalized 
flow c{l) exhibits the correct behavior near both the special 
and ordinary fixed points and in particular yields a canonical 
scaling of c with cross-over exponent (pord ~ -^v near the or- 
dinary transition. 
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I. INTRODUCTION 

Since the first renormalization group (RG) analysis of 
surface critical behavior by Lubensky and Rubin a 
number of subsequent advances ||^-|^ in the technique 
have enabled computation of exponents to O(e^), critical 
amplitudes, and various cross-over functions [|§-|]. In 
particular, Diehl and Dietrich have systematically devel- 
oped a formalism whereby the power and elegance of the 
field theoretic method has been fully exploited. 

For the case of an 0(1) system confined to the half- 
space z > it suffices to consider the reduced Hamiltio- 
nian 
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(1.1) 



the presence of the bounding surface manifesting itself 
as an additional surface interaction The parameter 
c takes account of the local enhancement of the reduced 
temperature in the vicinity of the surface. At lowest or- 
der the surface term results in the boundary condition 
(/)'(0) = c</)(0) and thus 1/c corresponds to the distance 
over which the order parameter falls to zero when ex- 
trapolated away from the surface. For c > the surface 
orders with the bulk while for c < there is an enhanced 
tendency to order at the surface. The "special" transi- 
tion with c = divides these two regimes whereas the 
"ordinary" transition with c = oo corresponds to a state 
where ordering at the surface is completely supressed. 

An issue of considerable interest is the manner in which 
various quantities behave close to the ordinary transition 
and to this end an expansion in the hare extrapolation 
length 1/c has been developed ||ll|.|lj] . Among the results 
is the finding that at the ordinary point energy related 



quantities involving (/)^ averages exhibit behavior charac- 
terized by relations involving hulk exponents [p^ . This 
in turn is a direct consequence of a vanishing anomalous 
exponent r\c associated with the extrapolation length. 

The canonical scaling of c(Z) near the ordinary point is 
interesting in light of the fact that all analyses address- 
ing the cross-over behavior from the special to ordinary 
point [0,^1^ have utilized a linear RG equation which in 
dimensional regularization reads 



dc 
dl 



(1 + r]c)c 



(1.2) 



with e corresponding to the block spin size. For finite 



c{l) Eq. (1.2) results from a straightforward application 
of the field theoretic method to a bulk system with a pla- 
nar bounding surface 1^,^. Clearly equation (1.2) yields 
a flow which is independent of the proximity to the spe- 
cial transition and thus does not display the expected 
classical behavior at large c{l). It is readily verified, how- 
ever, that this disparity in scaling behavior from that in- 
ferred from the 1/c expansion is compensated for by the 
crossover functions exhibiting logarithmic singularities at 
large c{l). When exponentiated, these singularities lead 
to powers of c{l) that in effect undo the incorrect large c 
behavior of Eq. ( |l.2| ) and in turn lead to the appropriate 
exponents at the ordinary point. 

With the above in mind a question of immediate inter- 
est is to what extent it is possible to deduce a flow for c(/) 
that correctly interpolates between the special and ordi- 
nary points. Constructing such a flow is not immediately 
obvious since for finite c the linear RG Eq. (1.2) results 



from the standard program of renormalizing all relevant 
surface operators. It has been pointed out ||ll|, however, 
that near the ordinary point additional care must be used 
in categorizing operators as relevant or irrelevant in the 
RG sense. In particular, in the context of the 1/c expan- 
sion it happens that insertions of the formally irrelevant 
interaction {dn<p)'^ must be considered ||ll[]. 

In the following we outline an approach based on the 
physical notion of an extrapolation length that gives rise 
to a non-linear RG equation exhibiting the correct be- 
havior both at the special and ordinary fixed points. As 
in the case of the 1/c expansion the operator (9„0)^ is 
found to play a important role in the analysis. In addition 
to yielding an RG flow with the sought-after behavior at 
both fixed points, the method further elucidates the con- 
nection between the extrapolation length and parameter 
c. This insight is of interest in its own right since conven- 
tional wisdom holds that the connection loses meaning 
beyond mean-field theory Q . Making use of an analysis 
based on phase shifts, we will demonstrate that there is a 
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means of extending the notion of an extrapolation length 
that remains intact when fluctuations are taken into ac- 
count. This approach may have the potential for further 
development since the use of phase shifts in many-body 
systems is a concept often incorporated into perturbative 
analyses of arbitrary order . 



II. SCATTERING FORMULATION 

O ur a nalysis begins with the reduced Hamiltonian of 
Eq. (1.1) with the undestanding that all volume integra- 
tions arc to be taken over the half-space z > 0. Within 
mean fi eld t heory, straightforward variation of the Hamil- 
tonian (LI) gives rise to the boundary condition 



cm = 0'(o) 



(2.1) 



Beyond mean field theory there is an approach in which 
the connection between c and an extrapolation length is 
still apparent. The key observation is that the oscillatory 
nature of the modes leads to the extrapolation length 
manifesting itself as a phase shift. In particular, taking 
the free Hamiltonian TLq to be 
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the modes which diagonalize Tio satisfy the boundary 
condition Eq. {2A) and are of the form 
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where a{z) = S(z) — S(oo), and t = r + I](oo) is the 
suitably shifted bulk reduced temperature. Since we will 
ultimately be implementing the R G, l et us assume that 
cr ~ e which then permits Eq. (2.6) to be solved us- 
ing standard perturbation theory. In the present circum- 
stance in which results appropriate to one-loop order are 
sought, first order perturbation theory suffices. We are 
led to a modified scattering amplitude 

fr = f-^l dz[cbl{z)?aiz) (2.7) 

The corresponding fluctuation corrected Cr is obtained 
by considering the small k limit of Eq. ( |2.5| ). Noting 
that / ~ 1 -I- 2ik/c it follows that 
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— = - dz{z + l/cf(j{z) 

Cr C J 



(2.8) 



which is valid to 0{e). In the event that higher order 
corrections are desired it is necessary to ta ke a ccount of 
additional perturbative corrections to Eq. (2.6). 



III. RENORMALIZATION GROUP 

The above results can now be used to determine the 
rcnormalization group equation for c(Z). We proceed with 
a momentum-shell approach. To this end, we note that 
translational invariance in directions parallel to the sur- 
face allows one to write 



(x)=^0,(z)e''i-> 



(3.1) 



where the scattering amplitude fk and phase shifts are 
given by 



fk — 



c + ik 
c — ik 



tanSk — k/c 



(2.4) 



At this level of approximation it is apparent that the 
presence of the surface interaction serves as an effec- 
tive scattering potential characterized by phase shifts Sk- 
Conversely, given phase shifts dk the associated inverse 
extrapolation length satisfies 
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The inclusion of fluctuations will alter the effective sur- 
face potential. To ascertain how fluctuations influence 
the extrapolation length one can appeal to the manner 
i n w hich the phase shifts are modifled and then use Eq. 

We now carry out this program by addressing the low- 
est order, one-loop corrections. For a given self energy E 
the modes satisfy 



dz^ 



+ it + q^ 



Eh 



-a{z) 



(2.6) 



Integrating out all modes with parallel momentum in the 
shell e~^' < q <1 and using the result that the averages 
obey 



-k\z — z'\ 



a e 



-k(z+z') 



with 



(3.2) 



(3.3) 



and = +t, one finds for the subtracted self energy 



a{z) ^ - 



uK, 



4ki 



aic 



(3.4) 



where the 1-subscript refers to all quantities being eval- 
uated at g = 1. Rescaling lengths so that c 
arrives at the RG equation 



C^'c, one 



dc u*Kd- 



C — Kl 



2kj (c + ki) 



(3.5) 



where u has now been set to its fixed point value 
u*Kd-i = 8e/3. For comparison we note that the cor- 
responding equation resulting from a standard momen- 
tum shell approach in which only the interactions 
d>^,(hdn(b are considered reads 
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dc u*Kd- 



(3.6) 



Equation (refeqircl) is the har d cu t-off version of the di- 
mensionally regularized result ( |l.2| ). This latter equation 
being linear in c implies a flow 



c{l) 



(3.7) 



that is independent of the proxi mity to s pecial transition. 
Deviations between equations (3.6, 3.5) begin to appear 
when c(l ) ^ 1. The third non-linear term appearing in 
Eq. (3.5) for finite c is ultraviolet convergent and corre- 
sponds to the inclusion of contributions from the formally 
irrelevant {dn4>)'^ vertex. However, if this last term is ex- 
panded in 1/c, it is clear that corrections to the shift and 
exponent of c occur, and that successive terms become 
increasingly ultraviolet divergent. 

Equation (^) leads to some interesting results, which 
we now address. For the sake of illustration assume that 
the system is close enough to criticality so that cross- 
over to the ordinary point has already occurred while 
r(l) <IC 1. In this case Eq. (|3.5|) reduces to 



dc 
dl 
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c- 1 



(c - 1) 
2 {c+1) 



(3.8) 



Although it is possible to solve Eq. ( |3.8| ) exactly, only 
results accurate to 0(e) will be considered. There are 



several ways to go about solving Eq. (3.J), one of which 
proceeds by iteratively solving the differential equation 
to 0{e). This leads to the explicit solution 



c{l) = bil) 



u*Kd-i 



1 
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b{l)\n{l + b{l)) 



(3.9) 



6(0 = 6(0)6(1+"'=)' 



with 5(0) = c(0) + u*Kd-i/8, and rj^ = -e/3. Another 
method approximates the r oot s to the resulting cubic on 
the right hand side of Eq. (3.S) and leads directly to the 
implicit form 



C{1) - IJc 
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f c{l) H 


-Vc' 




■2/77c H 


-c(0)1 




[c(0) - 77c J 
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-c(0)- 








Kc(Oj 


= e 



(3.10) 

wh ich c an also be shown to follow from exponentiation 
of ( 3.10D . Inspection of the above results reveals that for 
c{l) <^ 1 the flow is characterized by r/c 7^ while close 
to the ordinary point c{l) ~ e' thus implying a vanishing 
r]c- Stated differently, the cross-over exponent for the 
extrapolation length X{1) = l/c{l) at the ordinar y po int 
is (fiord = —V. Another interesting feature of Eq. (3.10) is 
that it yields an ordinary fixed point of order c{oo) ^ 1/e. 

Past analyses, employing the momentum-shell tech- 
nique to surface related phenomena, have encountered 
various technical difficulties Vim- We, therefore, first 



consider how this method leads to the standard linear 



equation (1.2) before attempting an alternate derivation 
of the modified RG equation ( |3.5| ) . As degrees of freedom 
are integrated out, additional interactions are generated. 
This is accommodated by taking the surface interaction 
to be of the form 



V{z) 



(3.11) 



with S^"^^ (z) referring to a m}^ derivative. For given V{z) 
the coefficients Vm are determined by 
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z"'V{z)dz 



(3.12) 



Consider the one-loop contribution, which results in the 
surface interaction V{z) — a{z) given by Eq. (3.4). Af- 
ter rescaling the surface spins by a factor e'^'*^^" 



mJ/2 



one 



finds that the coefficients 
tions: 



satisfy the recursion rela- 



dVm 

dl 



(1 



uKd^i (-)'"ai 
2 (2ki)™+2 



(3.13) 



The vertex involving (fPS'{z), or equivalently 5{z)(f>dn4) 
results from the boundary term associated with (V(^)^. 
Analogous to what is done in bulk phenomena the factor 
771 is chosen so that vi — 1/2 remains fixed. This leads 
to the result 



u*K, 



d-l 



-ai 



(3.14) 



When this value for 771 is inserted into Eq. ( |3.13D for 
Wo, one ends up with the linear RG equation ( p_.2| ). It 
is interesting that the non-linearity associated wit h the 
factor ai is entirely cancelled. Inspection of Eq. ( 3.13| ) 
governing Vm reveals that all interactions with m > 2 are 
irrelevant. However, in the context of calculating various 
scaling functions, these interactions with m > 2 must, in 
fact, be considered to account for all 0(e) contributions 

It is possible under certain circumstances to interpret 
the contribution to c from rji as feeding in from the vi 
vertex. This becomes evident upon considering the con- 
tribution each surface term makes when inserted into a 
propagator with legs off the surface. Recall that the m^^ 
vertex involves a factor (5(™-'(z), which leads, after an 
integration by parts, to an interaction 5(z)9™[0^]. The 
boundary condition then effectively relates this to a term 
proportional to S{z)(l)'^ and leads to a contribution from 
v„i feeding into the recursion for c. In such case it is 
possible that the higher order surface interactions will 
infiuence the behavior of c. 

To see how the above reasoning leads to the modified 
RG Eq. (3.5) assume that all two point interactions ulti- 
mately will modify a propagator with legs off the surface 
and that these legs each carry a transverse momentum q 
with associated factor kq ~ \/ q'^ + r{l). The assumption 
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that the legs are off the surface effectively leads to the 
replacements: 



(3.15) 
(3.16) 



in all surface interactions. Note that the 8 function sin- 
gularity associated with two or more derivatives makes 
no contribution because of the fact that the legs are off 
the surface. When momenta in a thin shell are integrated 
out, each Vm receives a contribution 



td (-)"ai 



(3.17) 



For the moment, we will ignore any contribution from 
an anomalous surface spin rescaling factor r]\ . Integrat- 
ing by parts and invoking the correspondence (3.16), the 
v\ interaction leads to a term Awi5(z)c0^. Similarly, 
the vi vertex involves (9^[(/)^] and thus leads to a term 



2At)2(5(z)(c-^ 



It follows that the total effective 



contribution from the vx , V2 interactions to is 
Auo = Awic -t- 2Au2(c^ + k§) 



(3.18) 



Rescaling spins and lengths, using Eq. (3.17), and for 
the moment ignoring the last term involving kq, one ar- 
rives at the modified RG Eq. (3.5). Alternatively, iden- 
tifying Aum with the moments of V{z) using Eq. (3.12) 



one recovers Eq. (2^) derived from the scattering theory 
approach. Though r/i was ignored, the final result is the 
same when anomalous spin rescaling is included. If spins 
are rescaled so that v\ remains fixed, while there is no 
contribution to vq in the form of Awi, there is a contri- 
bution from 771 which, because of Eq. (3.14), yields an 
identical result. 

The above analysis arbitrarily neglected the contribu- 
tions from the U2 vertex in addition to all interactions 
with m > 2. To determine under what circumstances 
this is justified note that the vertex w,„ makes a contri- 
bution to Awo of order 



IV. CONCLUDING REMARKS 

We have presented a method for identifying an effec- 
tive surface enhancement c(Z) which utilizes the scatter- 
ing phase shifts of the localized part of the self energy 
(t(z). The resulting cross-over behavior in c(/) is found to 
arise from the inclusion of contributions of various higher 
order surface interactions, in particular (9„(/))^. It is in- 
teresting that the relatively simple connection involving 
phase shifts implicitly includes such higher order correc- 
tions. Furthermore, the method appeals to characteris- 
tics of the entire (smooth) surface interaction rather that 
its constituent localized (delta function) pieces and thus 
may lead to further insights into surface phenomena. In- 
deed, though often convenient, the use of hyper-localized 
surface distributions is somewhat unphysical and occa- 
sionally leads to pathological quantities requiring special 
limiting procedures and interpretations. 

The task of determining a scaling field with the cor- 
rect scaling behavior at both fixed points is important in 
its own right. We have performed a preliminary analysis 
of the scaling functions for the surface susceptibility and 
surface free energy and find, as ex pecte d, that the use 
of a modified flow similar to Eq. ( 3.1C ) eliminates the 
logarithmic singularities otherwise found in these quan- 
tities [^,1). This in turn suggests that the logarithmic 
singularities in these two quantities are due entirely to 
the cross-over in ciV). 

Within his calculation of the local susceptibility Gold- 
schmidt also addressed the exponentiation of the log- 
arithmic singularities appearing in this quantity. In this 
particular scaling function, however, the introduction of 
our modified flow is not sufhcient to eliminate the singu- 
larity. We have also verified this is also the case for the 
layer susceptibility This is to be expected, however, 
since both these quantities involve at least one external 
point on the surface. 



[(c + ACo)' + (-)"(c-«;o)' 



(3.19) 



and becomes increasingly negligible for kq <IC ki 1. 
This latter condition is satisfied sufficiently close to the 
critical point when leg momenta q<^ 1. For the cur- 
rent situation of interest here this condition is well satis- 
fied. However, in view of this assumption, our derivation 
strictly applies only to the RG Equation (3.8), in which 
r{V) was neglected. 

Generally, when r[V) is not negligible it is possible to 
sum the higher order corrections that were neglected in 
the above analysis. However, the resulti ng e quation dif- 
fers from that found using phase shifts (pT^). The dif- 
ferences arise from the two methods reflecting different 
conventions on the finite part of c{V). This is, of course, 
compensated for by making a correspondingly different 
subtraction, depending on which flow is used. 
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